ELSEVIER 


Available  online  at  www.sciencedirect.com 

ScienceDirect 

Journal  of  Power  Sources  178  (2008)  537-546 

www.elsevier.com /locate /jpowsour 


Percolative  model  of  proton  conductivity  of 
Nation®  membranes 

Paola  Costamagna*,  Simone  Grosso,  Renzo  Di  Felice 

DICheP,  Department  of  Chemical  and  Process  Engineering  ‘G.B.  Bonino’,  University  of  Genoa, 
Via  Opera  Pia  15,  16145  Genoa,  Italy 

Received  18  July  2007;  received  in  revised  form  23  September  2007;  accepted  24  September  2007 
Available  online  29  September  2007 


Abstract 

A  model  is  proposed  for  the  simulation  of  Nation®  proton  conductivity,  where  it  is  assumed  that  proton  conduction  occurs  only  in  the  water  present 
in  the  membrane  pores.  Water  is  considered  to  be  present  in  the  pores  due  to  two  different  phenomena:  adsorption  and  capillary  condensation. 
In  the  latter  case,  the  pore  is  flooded  and  proton  conduction  occurs  throughout  the  whole  pore  section.  The  conditions  under  which  capillary 
condensation  occurs  are  simulated  in  the  model  through  the  Kelvin-Cohan  equation  for  condensation.  The  Kelvin-Cohan  equation  is  a  function 
of  RH,  temperature  and  the  pore  radius;  the  larger  the  pore,  the  higher  the  RH  for  which  capillary  condensation  takes  place.  If  the  conditions  for 
capillary  condensation  are  not  satisfied,  then  water  is  present  in  the  pore  due  to  adsorption  under  the  form  of  a  water  layer  which  covers  the  pore 
walls  and  provides  a  path  for  proton  conduction.  In  this  case,  the  modified  Brunauer-Emmet-Teller  (BET)  equation  has  been  used  in  the  model  to 
simulate  the  thickness  of  the  water  layer. 

In  both  cases  of  capillary  condensation  and  adsorption,  the  conductance  g  of  a  pore  has  then  been  calculated  through  the  formula  g  =  kSII,  where 
k  is  the  proton  conductivity  of  water,  S  the  cross-section  of  the  pore  volume  which  is  occupied  by  water,  and  /  is  the  pore  length.  Pores  of  different 
size  are  present  in  the  membrane  (data  of  pore  size  distribution  have  been  extracted  from  the  literature);  connectivity  of  the  water  layers  present  in 
the  different  pores  is  necessary  in  order  to  achieve  a  continuous  path  of  proton  conduction  through  the  membrane,  which  is  a  percolation  problem. 
To  this  end,  the  structure  of  the  membrane  pores  has  been  simulated  in  the  model  through  the  effective  medium  approximation  (EMA). 

The  simulation  results  of  proton  conductivity  of  the  membrane  show  good  agreement  with  literature  experimental  data,  even  when  varying  the 
RH  operating  conditions. 

©  2007  Elsevier  B.V.  All  rights  reserved. 

Keywords:  Transport  phenomena;  Modelling;  Nafion®;  Percolation  theory;  Porous  materials;  Fuel  cells 


1.  Introduction 

Nafion®  is  a  perfluorosulphonic  acid  membrane,  originally 
developed  for  application  in  the  chloro-alkali  process,  and 
now  applied  in  fuel  cells,  namely  proton  exchange  membrane 
fuel  cells  (PEMFCs)  and  direct  methanol  fuel  cells  (DMFCs). 
Nafion®  consists  of  a  perfluorinated  polymeric  backbone  (poly- 
tetrafluoroethylene,  PTFE)  and  side  chains  with  SO3H  groups. 
Several  studies  have  been  made  in  order  to  understand  the 
nanoscaled  structure  of  Nafion®,  since  the  early  work  of  Gierke 
et  al.  [1],  who  proposed  a  description,  still  generally  accepted. 
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where  the  polymeric  membrane  is  pictured  in  terms  of  an 
inverted  micellar  structure  in  which  the  ion  exchange  sites 
separate  from  the  fluorocarbon  backbone  thus  forming  spher¬ 
ical  clusters  (pores),  which  are  connected  by  channels,  with 
the  SO3H  groups  distributed  on  the  inner  surface  and  pointing 
inwards  (Fig.  1). 

When  the  membrane  is  humidified,  water  penetrates  the 
pores,  and  then,  due  to  their  acidic  nature,  the  sulphonic  groups 
dissociate  into  SO3-  groups,  which  remain  attached  to  the  back¬ 
bone  chains,  and  protons,  which  float  in  the  aqueous  phase,  and 
are  free  to  move  under  the  action  of  an  electrical  field.  Actually, 
it  is  ascertained  that  proton  conduction  is  not  occurring  through 
movement  of  the  same  proton,  but  instead  by  exchange  of  a 
free  proton  (i.e.  proton  hop)  between  one  H+-«H20  cluster  and 
another.  The  proton  hop  is  preceded  by  suitable  configurational 
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Nomenclature 

c  parameter  in  Eq.  (5) 

C  concentration  (mol  m- 3 ) 

EW  equivalent  weight  (kg  kmol- 1 ) 

g  conductance  (£2-1) 

gM  effective  network  conductance  (£2-1) 

k  parameter  in  Eq.  (6) 

l  pore  length  (m) 

nj  coordination  number 

p  pressure  (Pa) 

p°  saturation  pressure  (Pa) 

RH  relative  humidity 

r  pore  radius  (m) 

t  layer  thickness  (m) 

v  pore  volume-fraction 

Vl  molar  volume  of  bulk  liquid  (m3  mol- 1 ) 
x  specific  mass  of  adsorbed  gas  (kg  kg-1) 

xm  specific  mass  of  adsorbed  gas  in  a  monolayer 
(kg  kg-1) 

Greek  symbols 

y  surface  tension  (N  m- 1 ) 

e  void  degree 

k  proton  conductivity  of  water  (fi- 1  cm-1) 

X  moles  of  water  per  mole  of  sulphonic  group  in 

membrane 

Pmemb  density  of  dry  membrane  (kg  m-3) 
a  molecular  diameter  (m) 

(p,  i/  probability  density  distribution 

limit  molar  conductivity  (Q- 1  m2  mol- 1 ) 


polarization  fluctuations  in  water.  Subjects  of  investigation  [2,3] 
have  been  (i)  the  size  of  the  typical  H+-hHtO  clusters  (H30+, 
Hj02+,  H<)04+  or  a  mixture  of  them)  and  (ii)  the  different  types 
of  proton  conductivity  occurring  in  the  middle  and  on  the  surface 
of  a  water-filled  pore.  Indeed,  in  the  middle  of  the  water-filled 
pore,  proton  conductivity  resembles  proton  conductivity  in  pure 
free  water,  while  in  proximity  of  the  surface  of  the  pore  itself 
the  presence  of  the  SO3-  groups  modifies  the  proton  conduc- 


Fig.  1.  Structure  of  pores  in  Nation®  membranes.  (•)  Represents  the  — SO3H 
group. 


tion  mechanism.  In  this  respect,  Okada  et  al.  [4]  demonstrated 
that  some  of  the  water  molecules  in  the  membrane  are  strongly 
associated  with  the  SO3-  sites,  forming  the  so-called  primary 
hydration  layer  present  on  the  pore  surface;  the  ability  of  these 
water  molecules  to  form  hE-nFEO  clusters  and  provide  the  con¬ 
figurational  flexibility  necessary  in  order  to  contribute  to  proton 
conduction  is  highly  compromised,  resulting  in  very  low  proton 
conductivity  [5], 

Of  course,  in  order  to  allow  proton  conduction,  not  only  must 
water  be  present  in  a  pore,  but  also  the  aqueous  phase  must  per¬ 
colate  through  the  membrane,  forming  a  continuous  path  for 
proton  conduction.  Two  aspects  contribute  to  the  formation  of 
a  connected  network  of  water-filled  pores  inside  the  membrane, 
i.e.  the  amount  of  water  present  in  the  membrane  and  the  struc¬ 
ture  and  pore  size  distribution  of  the  pore  network.  Both  the 
previous  aspects  have  been  subject  of  experimental  and  theoret¬ 
ical  studies.  In  particular,  concerning  the  pore  size  distribution 
of  the  pore  network,  Divisek  et  al.  [6]  have  investigated  in  depth 
the  pore  size  distribution  of  Nation®  1 17  membranes,  by  apply¬ 
ing  a  method  of  standard  porosimetry  (SPM).  According  to  their 
results,  a  wide  spectrum  of  pore  diameters  was  detected,  ranging 
from  1  to  100  nm,  with  the  average  value  in  the  order  of  some 
nm.  Concerning  the  membrane  water  content,  several  experi¬ 
mental  studies  of  the  isotherms  of  water  adsorption  in  Nafion® 
are  reported  in  the  literature.  Indeed,  the  water  uptake  of  Nafion® 
1 17  membranes  has  been  measured  by  several  scientists  [7-1 1]; 
in  all  these  experiments,  water  adsorption  has  been  recorded,  by 
suspending  dry  samples  of  Nafion®  in  an  inert  environment  with 
controlled  relative  humidity  RH  (defined  as  the  ratio  between  the 
partial  pressure  of  water  vapour  p  and  the  saturation  pressure  p° 
at  the  operating  temperature);  water  uptake  is  then  measured  as 
the  difference  in  weight  between  the  dry  state  and  the  state  of 
equilibrium  with  water  vapour.  Unfortunately,  no  studies  have 
examined  both  adsorption  and  desorption  with  the  same  sam¬ 
ples.  Indeed,  only  one  study  has  been  found  on  desorption  [6]. 
The  desorption  data  have  been  obtained  after  soaking  the  mem¬ 
brane  in  liquid  water;  during  the  experiments,  the  membrane 
has  been  kept  in  a  gaseous  environment  with  decreasing  values 
of  RH,  starting  from  RH=  1;  equilibrium  was  obtained  before 
recording  any  data  point. 

The  experimental  data  obtained  by  the  different  authors  are 
reported  in  Fig.  2.  Since  the  different  membranes  were  not  in  the 
same  state  before  the  experiments  were  performed,  it  is  difficult 
to  discuss  the  differences  among  the  different  results  in  depth. 
However,  it  is  clear  that  desorption  and  adsorption  isotherms 
show  hysteretical  behaviour.  The  discrepancy  in  the  values  of 
water  impregnation  in  Nafion®  membranes  at  RH  =  1  obtained 
by  soaking  the  membrane  in  liquid  water  or  by  water  adsorption 
from  water  vapour  at  RH  =  1  is  well  known  as  the  Schroeder 
paradox  [9,11,12],  The  hysteretical  behaviour  of  Nafion®  with 
hydration/dehydration  cycles  is  also  mentioned  in  other  works 
[5,13],  Generally  speaking,  hysteresis  effects  can  be  explained 
by  the  existence  of  phase  transitions  [14—16],  or  by  the  visco¬ 
elastic  properties  of  Nafion®  membranes  [13],  or  by  the  classic 
vapour-liquid  phase  transitions,  that  is  capillary  condensation  in 
pores  with  diameters  ranging  from  1  to  10  nm,  which  is  exactly 
the  pore  size  of  Nafion®  as  measured  experimentally  by  Divisek 
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Fig.  2.  Experimental  literature  data  of  water  adsorption  and  desoption  isotherms 
in  Nation®  1 17  at  20-30  °C:  (A)  adsorption  at  25  °C  [8];  (■)  adsorption  at  25  °C 
[7];  (+)  adsorption  at  30  °C  [9];  (♦)  desorption  at  20  °C  [6]. 

et  al.  [6].  In  this  paper,  this  latter  interpretation  is  followed  and 
the  visco-elastic  properties  of  Nation®  membranes  are  not  taken 
into  account.  Thus,  we  assume  that,  when  starting  from  a  dry 
membrane  and  by  increasing  RH,  membrane  pore-filling  is  due 
to  water  adsorption  on  the  pore  walls  coupled  to  capillary  con¬ 
densation;  for  each  value  of  RH  some  of  the  pores  (the  smallest 
ones)  are  completely  filled  with  water,  while  the  largest  pores 
contain  only  a  thin  layer  of  water  adsorbed  on  their  pore  surface. 
Thus,  the  conductance  of  each  pore  is  thus  very  different  depend¬ 
ing  on  its  radius,  and  the  overall  conductivity  of  the  membrane 
depends  on  the  percolation  through  the  pore  network. 

In  this  paper,  a  model  of  Nation®  membranes  is  proposed, 
which,  starting  from  the  pore  size  distribution  proposed  by 
Divisek  et  al.  [6],  evaluates  the  occurrence  of  capillary  con¬ 
densation  in  the  pores  of  the  membrane,  calculates  their  proton 
conductance,  and  finally  the  overall  proton  conductivity  of  the 
membrane  taking  into  account  percolation  within  the  pore  net¬ 
work  through  the  effective  medium  approximation.  Each  of  these 
aspects  will  be  illustrated  in  more  detail  in  Section  2.  Section 
3  includes  a  comparison  between  model  results  and  literature 
experimental  data  of  water  adsorption  and  proton  conductivity  of 
Nation®  117  membranes  at  20  °C  and  at  different  values  of  RH. 
Indeed,  literature  studies  [17]  have  demonstrated  that  not  only 
the  water  uptake,  but  also  the  proton  conductivity  of  Nation® 
shows  a  dramatic  dependence  upon  the  operating  conditions  of 
relative  humidity  (RH),  with  several  orders  of  magnitude  of  dif¬ 
ference  between  the  values  obtained  at  10-20%  RH  and  at  100% 
RH  (where  maximum  conductivity  is  achieved).  This  is  a  signifi¬ 
cant  issue  for  fuel  cell  applications,  since  any  zone  of  the  fuel  cell 
where  the  relative  humidity  is  low  shows  a  dramatic  loss  of  per¬ 
formance.  However,  we  should  point  out  that  there  is  a  difference 
in  the  operating  conditions  of  a  Nation®  membrane  when  oper¬ 
ating  as  the  electrolyte  in  a  PEMFC,  and  the  Nation®  membrane 
simulated  in  this  paper.  When  a  Nation  membrane  is  embedded 
into  a  PEMFC  and  this  is  put  into  operation,  usually  the  level 
of  hydration  along  the  membrane  thickness  is  not  uniform.  The 
main  reason  is  that  the  RH  on  the  two  faces  of  the  membrane 
are  significantly  different  (in  particular,  the  RH  is  higher  at  the 
cathodic  side  due  to  the  water  released  by  the  electrochemical 


reaction,  and  due  to  the  protons  releasing  the  water  which  they 
have  dragged  through  the  membrane  due  to  the  electro-osmotic 
drag  phenomenon;  at  the  anodic  side,  on  the  contrary,  water  is 
continuously  dragged  away  by  the  protons  moving  towards  the 
cathodic  side).  Even  if  the  anodic  and  cathodic  feeding  flows 
into  the  fuel  cell  have  the  same  RH,  nevertheless,  these  feeding 
flows  are  usually  not  large  enough  to  allow  the  assumption  that 
the  water  released  at  the  cathodic  side  and  the  water  subtracted  at 
the  anodic  side  are  negligible  contributions.  As  a  consequence, 
the  value  of  RH  on  the  anodic  and  the  cathodic  face  of  the  mem¬ 
brane  is  different,  and  the  water  content  at  each  membrane/gas 
boundary  is  different  as  well.  This  causes  a  profile  of  membrane 
humidification  across  the  membrane  thickness,  which  in  turns  is 
the  driving  force  for  a  water  diffusion  flow  [  1 8] .  On  the  contrary, 
in  the  present  paper  we  have  considered  the  membrane  exposed 
to  the  operating  conditions  used  for  water  uptake  measurements 
[6-11]  and  for  the  proton  conductivity  measurements  [17],  i.e. 
the  membrane  is  suspended  in  an  environment  at  uniform  and 
controlled  temperature  and  relative  humidity.  In  our  model,  both 
sides  of  the  membrane  are  in  identical  conditions  and  at  the 
same  RH,  and  thus  the  membrane  humidification  (which  is  eval¬ 
uated  as  an  average  of  the  water  content  of  the  pores  of  different 
dimensions  present  in  the  membrane  itself)  is  the  same  at  each 
gas/membrane  interfaces  and  also  throughout  the  whole  mem¬ 
brane,  and  thus  no  water  diffusion  is  expected  to  take  place.  In 
addition,  in  PEMFCs  the  number  of  water  molecules  dragged 
by  the  protons  along  the  membrane  due  to  the  electro-osmotic 
drag  phenomenon  varies  along  the  membrane  thickness  since 
it  depends  on  the  local  value  of  the  membrane  hydration  [18]; 
the  number  of  water  molecules  dragged  by  the  protons  is  lower 
in  the  regions  where  membrane  hydration  is  poorer.  This  is  a 
further  mechanism  of  water  transport  through  the  membrane 
since  in  this  way,  for  example,  if  protons  move  from  a  region 
at  high  humidity  to  a  region  at  low  humidity,  they  transport 
water  from  the  previous  to  the  latter  region.  This  effect  is  not 
considered  in  the  present  model  since,  as  explained  above,  the 
hydration  level  of  the  membrane  is  uniform  through  the  whole 
membrane. 

Thus,  the  model  presented  in  this  paper  does  not  include  water 
diffusion  and  the  electro-osmotic  drag,  since,  in  the  operating 
conditions  considered  here,  the  first  is  not  expected  to  occur  and 
the  second  is  not  expected  to  have  an  effect.  For  the  purposes 
of  PEMFC  simulation,  the  present  model  can  be  coupled  to  the 
model  of  water  diffusion  and  electro-osmotic  drag  through  the 
membrane. 


2.  Model  formulation 

The  model,  starting  from  a  pore  size  distribution  of  the 
Nation®  membrane,  evaluates  the  occurrence  of  capillary  con¬ 
densation  in  the  pores  of  the  membrane,  calculates  their  proton 
conductance  and  finally  the  overall  proton  conductivity  of  the 
membrane  taking  into  account  percolation  within  the  pore  net¬ 
work  through  the  effective  medium  approximation.  The  model 
is  structured  as  depicted  in  Fig.  3.  In  the  following,  each  step  is 
described  in  detail. 
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2.1.  Pore  size  distribution 


The  pore  size  distribution  of  a  Nation®  117  membrane  has 
been  evaluated  experimentally  by  Divisek  et  al.  [6],  by  means 
of  a  method  of  standard  porosimetry.  In  this  method  (which  is 
considered  by  some  scientists  more  suitable  for  ridged  mate¬ 
rials,  such  as  the  gas  diffusion  layer,  than  for  membranes), 
they  were  able  to  obtain  curves  representing  the  integral  pore 
volume-fraction  distribution  v  (expressed  as  cm3  cm-13  of  over¬ 
all  membrane  volume),  in  terms  of  radii  for  Nation®  117 
membranes  at  20  and  80  °C;  the  data  at  20  °C  are  reported  in 
Fig.  4.  The  meaning  of  the  experimental  points  is  that  u;+i  —  u; 
is  the  overall  volume-fraction  of  all  the  pores  having  radius  rang¬ 
ing  between  r,  and  r;+i .  Under  the  assumption  that  all  the  pores 
have  the  same  length  (which  is  also  the  basic  assumption  behind 
the  Effective  Medium  Approximation  theory,  see  section  below), 
the  number  of  pores  having  a  certain  radius  r  is  proportional 
to  their  overall  volume  divided  by  r2;  thus,  the  experimental 
curves  of  v  can  be  used  to  obtain  the  pore  probability  density  by 


Fig.  4.  (■)  Probability  density  function  of  pore  radius  distribution  i/r(r)  for 
Nation®  1 17  at  20  °C,  derived  from  (•)  literature  experimental  data  of  integral 
pore  volume  distribution  v  at  20  °C  [6]. 


■  (cm3cm‘3) 
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applying  the  following  equation  to  the  experimental  points: 


(1) 


The  resulting  probability  density  function  i//(V)  is  then 
obtained  by  linearly  connecting  the  points  obtained  and  nor¬ 
malizing  the  area  to  1 .  The  results  obtained  for  the  probability 
density  function  evaluated  at  20  °C  are  reported  in  Fig.  4;  it  can 
be  noticed  that  it  is  very  narrow,  with  no  pores  below  1.05  nm 
and  above  about  10  nm.  Also,  it  displays  a  bimodal  behaviour 
with  two  peaks,  respectively  at  1.05  and  2.9  nm.  This  picture 
of  the  membrane  pores  is  in  quite  good  agreement  with  previ¬ 
ous  literature  results,  such  as  ones  reported  by  Gierke  et  al.  [1], 
which  proposed  a  widely  accepted  description  of  the  polymeric 
membrane  in  terms  of  an  inverted  micellar  structure  in  which 
the  ion  exchange  sites  separate  from  the  fluorocarbon  backbone 
thus  forming  clusters  (pores);  they  reported  that  when  the  mem¬ 
brane  is  in  the  dry  form,  an  average  cluster  has  a  radius  of  about 
1.8  nm. 


2.2.  Capillary  condensation  in  Nafion®  pores 

The  isotherm  reported  in  Fig.  5  is  referred  to  as  ‘Type  IV’ 
isotherm  [14],  It  is  characteristic  of  some  particular  types  of 
porous  solids  and  it  is  generally  accepted  that  it  appears  only  in 
solids  that  have  pore  diameters  ranging  from  tens  to  hundreds 
of  Angstrom  units.  Indeed,  after  the  physical  adsorption  process 
(branch  AB  in  Fig.  5),  that  is  identical  for  both  adsorption  and 
desorption  branches,  a  differentiation  occurs  between  adsorp¬ 
tion  and  desorption  branches.  In  the  adsorption  branch,  capillary 
condensation  occurs  (branch  BCE  in  Fig.  5)  due  to  instability  of 
the  film  interface  due  to  the  wall  curvature  arising  for  an  overall 
adsorption  layer  of  thickness  t  in  pores  of  radius  r.  This  causes 
water  condensation  in  liquid-like  conditions  leading  to  complete 
obstruction  of  the  pores  of  radius  r.  This  phenomenon  is  referred 
to  as  capillary  condensation;  the  classical  theoretical  approach 
used  to  describe  it  is  the  Kelvin-Cohan  equation  for  adsorption 

[19]: 

flrin(RH)  =  — —  (2) 

rc  -  t 


Fig.  5.  Type  IV  isotherm  of  porous  solids  (IUPAC  classification,  [14]). 


where  y  is  the  surface  tension,  Vl  the  molar  volume  of  the 
bulk  liquid  and  rc  is  the  pore  of  the  ‘critical’  radius:  pores 
with  radii  below  rc  are  completely  flooded  with  water.  Capillary 
condensation  is  characterized  by  hysteretical  behaviour,  since 
evaporation/desorption  follows  a  different  mechanism  (branch 
EDB  in  Fig.  5),  related  to  the  formation  of  a  meniscus  between 
the  condensed  fluid  and  the  vapour,  which  stretches  across  the 
pore;  in  this  case,  the  fluid  on  to  the  pore  walls  can  be  assimilated 
to  a  layer  of  thickness  t  and  follows  the  Kelvin-Cohan  equation 
for  desorption  [19]: 

f?rin(RH)  =  — t  (3) 

rc  -  t 

In  Eqs.  (2)  and  (3),  the  value  of  t  can  be  expressed  as  a  function 
of  RH  through  the  following  equation  [14]: 

,=(£)<’  (4) 
where  a  can  be  assumed  to  be  the  molecular  diameter  of  the 
adsorbate  (for  water,  a  %  3A),  and  xlxm  is  given  by  the  modified 
BET  equation  (see  Eqs.  (6)  and  (7)). 

Even  if  the  Kelvin-Cohan  equations  are  slightly  approx¬ 
imated  for  pores  smaller  than  20  nm  [14,20],  they  are  still 
currently  applied  in  a  number  of  applications,  and  for  example 
they  constitute  the  basis  for  the  Barrett-Joyner-Halenda  (BJH) 
method  for  the  pore  size  distribution  analysis,  which  is  imple¬ 
mented  as  the  standard  technique  on  all  automated  commercial 
adsorption  instruments. 

2.3.  Water  adsorption  in  Nafion®  pores 


The  process  of  gas/vapour  adsorption  in  porous  solids  and 
membranes  has  been  extensively  studied  [14,20,21],  and  is  the 
basis  of  experimental  techniques  leading  to  the  estimation  of 
pore  size  distribution  and  specific  surface  area.  In  the  general 
case  of  a  non-polar  adsorbate  and  a  porous  solid  with  no  iono- 
genic  groups,  the  process  occurring  at  low  partial  pressures  of  the 
adsorbate  in  the  gas  phase  in  thermodynamic  equilibrium  with 
the  porous  solid  is  the  physical  adsorption.  Physical  adsorption 
is  defined  as  the  condensation  of  gases  on  free  surfaces  [14] 
and  is  considered  to  involve  as  a  first  step  in  the  formation  of 
a  monolayer,  i.e.  a  completely  filled  single  layer  of  adsorbate 
molecules  accommodated  on  the  accessible  internal  surface  of 
the  porous  solid.  By  increasing  the  partial  pressure  of  the  adsor¬ 
bate  in  the  gas  phase,  further  layers  start  to  build  up  on  to  the 
first  monolayer  (multilayer  adsorption).  The  equation  generally 
used  to  describe  multilayer  adsorption  in  porous  solids  is  the 
Brunauer-Emmet-Teller  (BET)  equation,  which  is  derived  from 
a  kinetic  approach  to  the  process  of  adsorption: 


*m  (1  -  P/P°)W  +  (c  -  1  )p/p°] 

where  x  is  the  specific  mass  of  the  adsorbed  gas,  xm  the  specific 
mass  of  the  adsorbed  gas  on  the  monolayer  and  c  is  a  parameter 
related  to  the  molar  enthalpy  of  adsorption,  c  and  xm  are  usually 
treated  as  fitting  parameters,  which  are  evaluated  by  fitting  the 
experimental  data  of  adsorption  through  Eq.  (5). 
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However,  when  the  adsorbate  is  a  polar  molecule  (as  water) 
and  ionogenic  groups  are  present  in  the  membrane,  the  process 
of  adsorption  is  primarily  due  to  solvation  of  the  membrane 
ionogenic  groups,  and  the  BET  equation  is  not  applicable  tout 
court.  As  reviewed  by  Freger  et  al.  [22],  in  this  case  several 
approaches  have  been  used  to  obtain  an  analytical  equation  for 
the  water  uptake  as  a  function  of  water  activity,  that  is  the  water 
sorption  isotherm.  Indeed,  the  BET  equation  has  been  applied, 
under  the  hypothesis  of  treating  the  phenomenon  of  ion  solva¬ 
tion  as  a  typical  adsorption  process,  but  it  has  been  demonstrated 
to  be  reliable  only  for  low  activities,  up  to  ~0.3.  Other  equations 
have  been  proposed  and  among  them  the  modified  BET  equa¬ 
tion  has  been  demonstrated  to  be  suitable  for  the  whole  range 
of  activities  spanning  from  0  to  1.  This  model  is  widely  applied 
and  commonly  accepted  for  the  interpretation  of  isotherms  of 
water  adsorption  in  hydrophilic  polymers  [22],  It  is  based  on  the 
consideration  that  the  activity  appearing  in  the  BET  equation 
must  be  multiplied  by  an  ‘activity  coefficient’  k,  which  accounts 
for  the  particular  type  of  water-water  interaction  in  the  vari¬ 
ous  solvation  layers.  Thus,  the  modified  BET  equation  reads  as 
follows: 

X  _  ck(p/pO) 


Table  1 

Value  of  model  parameters 


RH[~] 

Fig.  6.  ( — )  Interpolation  of  experimental  literature  data  of  water  adsorption  and 
desorption  isotherms  for  Nafion®  1 17  at  20-30  °C  in  the  region  of  superimposi¬ 
tion,  through  the  modified  BET  equation  (Eq.  (7));  (A)  adsorption  at  25  °C  [8]; 
(■)  adsorption  at  25  °C  [7];  (+)  adsorption  at  30  °C  [9];  (♦)  desorption  at  20  °C 
[6], 


Model  parameters  (at  20  °C) 

Value 

l 

21.5  nm 

Vl 

18  g  mol-1 

y 

73dyncm_1 

K 

2.38  S  cm-1 

a 

3  A 

V  (Eq.  (7)) 

2.3  mol(H20)/mol  (SO3-) 

c  (Eq.  (7)) 

70 

k  (Eq.  (7)) 

0.7 

xm  (1  -  kp/p°)[  1  +  (c  -  1  )kp/p°] 

where  the  number  of  fitting  parameters  is  three;  the  physical 
meaning  of  the  parameters  xm  and  c  is,  in  principle,  the  same  as 
that  of  the  BET  equation. 

Considering  the  experimental  points  reported  in  Fig.  2, 
the  branches  of  the  adsorption  and  desorption  curves  which 
superimpose  can  be  considered  as  the  branches  where  only 
adsorption/desorption  occurs,  and  thus  that  part  has  been 
interpolated  through  the  modified  BET  equation,  rewritten  as 
follows: 

A  =  ck(p/p°) 

km  (1  -  kp/p°)[  1  +  (c  -  1  )kp/p°] 


Fig.  6  reports  the  results  of  the  interpolation;  the  values  of  the 
parameters  ).m,  c,  k  derived  from  a  least  square  fitting  method 
are  reported  in  Table  1 . 

2.4.  Evaluation  of  proton  conductance  of  the  pores 

The  model  is  based  on  the  hypothesis  that  proton  conduc¬ 
tion  occurs  within  the  water  contained  in  each  pore  of  the 
Nafion®  membrane.  This  approach  is  proposed  and  followed  by 
a  number  of  researchers  [5]  (more  comprehensive  approaches 
are  reviewed  in  [23]).  Thus,  a  simplified  evaluation  of  the  proton 
conductivity  of  the  water  contained  in  the  pores  of  the  membrane 
can  be  made,  based  upon  the  evaluation  of  the  concentration  of 
protons  into  the  water  contained  in  the  membrane  pores,  which, 
in  turn,  is  obtained  on  the  basis  of  the  equivalent  weight  of  the 
membrane  (EW=  1 100  g  mol-1),  the  density  of  the  dry  mem¬ 
brane  (pmemb  =  2.24  g  cm-3)  and  of  the  overall  volume-fraction 
of  the  pores  of  the  membrane  (v  ~  0.4,  see  Fig.  4),  in  the  fol¬ 
lowing  way: 


C  = 


Pmemb 

u-EW 


(8) 


The  concentration  calculated  in  this  way  is  C  =  5.1  M.  Under 
these  conditions  of  high  concentration,  the  molar  conductivity 
can  be  assimilated  to  the  limit  molar  conductivity  of  protons  in 
water,  which  is  A ^  =  350  S  cm2  mol-1  at  ambient  temperature 
[24],  Thus  the  conductivity,  calculated  as: 

k  =  A°cC  (9) 

is  2.38  S  cm-1  at  20  °C.  This  value  is  then  used  to  calculate  the 
conductance  of  each  pore,  keeping  in  mind  that  for  every  value 
of  the  relative  humidity  the  water  content  of  each  pore  depends 
on  the  pore  radius.  For  pores  with  radius  smaller  than  the  crit¬ 
ical  radius  rc,  satisfying  the  Kelvin-Cohan  equation  (Eq.  (2)) 
water  condensation  in  liquid-like  conditions  leads  to  complete 
obstruction  of  the  pores.  Thus,  for  a  given  temperature  and  rel¬ 
ative  humidity,  all  the  pores  having  radius  smaller  than  rc  are 
completely  filled  with  water,  while  pores  having  a  radius  larger 
than  rc  only  have  their  walls  covered  by  a  water  layer,  as  depicted 
in  Fig.  7.  The  thickness  of  the  water  layer  in  pores  is  calculated 
from  Eq.  (4),  where  er  is  the  molecular  diameter  of  water. 

Under  these  conditions,  once  the  water  content  of  each  pore  is 
calculated,  its  proton  conductance  is  easily  evaluated  by  apply¬ 
ing  Ohm’s  law.  For  a  pore  completely  filled  with  water,  the 
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Fig.  7.  Schematisation  of  pores  with  different  degrees  of  water  filling:  (a)  surface 
adsorption  (r>  rc);  (b)  capillary  condensation  (r  <  rc).  Arrows  represent  protons. 


proton  conductance  g  is: 


(10) 


If  only  a  layer  is  present,  the  following  formula  is  applied: 


2jtrt 

Knr 


(11) 


Actually,  in  order  to  take  into  account  that  the  first  layer  of 
water  molecules  attached  to  the  membrane  pores  are  somehow 
bridged  to  the  SO3-  groups  and  thus  not  completely  free  to 
rotate  and  allow  proton  conduction,  water  layer  thickness  thinner 
than  3.8  A  have  been  considered  not  to  contribute  to  the  proton 
conduction  mechanism. 

Starting  from  the  previously  mentioned  distribution  of  pore 
sizes,  the  probability  density  function  of  the  conductance  within 
the  membrane  is  also  evaluated  (applying  the  rules  for  the  func¬ 
tions  of  random  variables): 


, ,  ,  f(.gl/xK)  t((gl/2mct)  +  (f/2)) 

«s)=\2vsm+ — 12^74 —  <12) 

When  RH  is  low,  the  water  film  adsorbed  on  the  surface  walls 
of  the  pores  is  thin  (Fig.  7a),  the  contribution  of  these  pores  to 
the  overall  proton  conduction  is  small,  and  the  proton  conduc¬ 
tivity  of  the  membrane  is  based  only  on  the  transport  occurring 
in  the  flooded  pores  (Fig.  7b).  Under  these  conditions,  connec¬ 
tivity  among  the  flooded  pores  in  the  membrane  is  a  key  aspect. 
Indeed,  when  the  concentration  of  the  flooded  pores  is  low,  they 
are  either  isolated  or  grouped  in  small  isolated  clusters.  When 


RH  is  increased,  then  concentration  of  flooded  pores  increases 
and  larger  clusters  are  formed,  i.e.  the  mean  size  of  clusters 
increases  monotonically.  As  the  concentration  of  flooded  pores 
approaches  the  percolation  threshold  from  below,  the  larger  clus¬ 
ters  begin  to  merge,  creating  a  few  extremely  large  clusters,  so 
that  at  the  percolation  threshold  the  mean  cluster  size  diverge  to 
infinity,  which  means  that  for  any  size  of  the  system  (membrane) 
under  consideration,  there  is  a  complete  path  of  adjacent  flooded 
pores  crossing  the  membrane,  and  thus  macroscopic  proton  flow 
becomes  possible. 


2.5.  Simulation  of  the  pore  network  through  the  effective 
medium  approximation  (EMA) 


The  EMA  has  been  used  in  order  to  evaluate  the  conductiv¬ 
ity  of  the  disordered  pore  structure  within  the  membrane.  This 
method  was  first  proposed  to  evaluate  the  electrical  permittivity 
of  a  binary  random  mixture  of  continuous  phases  [25],  and  later 
on  it  has  been  applied  to  the  evaluation  of  the  conductivity  of 
random  conductor  networks  [26,27].  In  this  approach,  the  dis¬ 
ordered  medium  is  replaced  by  a  uniform  system  (Fig.  8)  whose 
properties  are  described  by  an  analytical  equation  [26,27]: 


(gM  ~  g)  0(g) 
g  +  (nT/ 2-  1)  gM 


dg  =  0 


(13) 


Eq.  (13)  allows  the  evaluation  of  the  effective  conductance  gM 
of  the  uniform  system  consisting  of  nodes,  each  connected  to 
its  neighbours  by  a  number  of  conducting  elements  (=2,3,4, 
...  00),  having  conductance  ga<g<gb;  g  varies  randomly 
according  to  a  probability  density  0(g).  Finally,  the  effective 
conductance  gM  evaluated  for  the  uniform  system  is  also  an 
estimate  for  the  property  of  the  starting  disordered  medium.  The 
EMA  has  been  demonstrated  to  be  very  accurate  if  the  system 
is  far  from  the  percolation  threshold,  regardless  of  the  expres¬ 
sion  used  for  0(g)  (however,  the  EMA  predictions  may  be  less 
accurate  close  to  the  percolation  thresholds)  [28], 


forming  to  the  EMA  assumptions. 
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r(nm) 

Fig.  9.  ( — )  Probability  density  function  of  pore  radius  distribution  data  for 
Nafion®  117  at  20  °C  (extracted  from  [6]).  ( — )  Calculation  of  the  relative 
humidity  RH  in  correspondence  of  which  a  certain  radius  becomes  critical  for 
capillary  condensation  at  20  °C. 


Fig.  10.  Pore  conductance  g  as  a  function  of  pore  radius  r 
20  °C  for  various  values  of  relative  humidity  RH. 


2.6.  Pore  length  estimation 

The  evaluation  of  the  pore  size  distribution  has  also  allowed 
estimation  of  the  typical  pore  length;  indeed,  conforming  to  the 
assumptions  made  for  the  EMA,  by  assimilating  the  pore  net¬ 
work  structure  to  a  regular  cubic  network,  the  following  equation 
applies: 


'2f(r)dr 


Considering  that  the  void  degree  of  the  membrane  is  about 
e  ~  0.4  (see  also  Fig.  4),  the  estimated  pore  length  is  about  20  nm. 


3.  Results  and  discussion 

The  model  is  independent  on  membrane  thickness,  and  thus 
the  results  presented  in  this  section  hold  for  Nafion®  membranes 
of  any  type  (111,  112,  115,  117,  etc.),  provided  that  their  pore 
size  distribution  is  the  same  as  that  used  in  this  paper  [6],  which 
was  obtained  specifically  from  Nafion®  117  membranes. 

Fig.  9  reports  the  probability  density  function  of  pore  radius 
distribution  for  a  Nafion®  117  membrane  at  20  °C;  the  same 
figure  also  shows  the  relative  humidity  RH  in  correspondence 
of  which  a  pore  having  a  certain  radius  becomes  completely 
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Fig.  11.  Model  results  of  probability  density  of  pore  proton  conductance  g  in  Nafion®  1 17  at  20  °C  at  different  values  of  relative  humidity  RH. 
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flooded  due  to  capillary  condensation  (according  to  Eq.  (2)). 
Fig.  9  shows  that  the  smallest  pores  of  the  membrane  (radius 
of  1.05  nm)  become  completely  flooded  with  water  at  a  rel¬ 
ative  humidity  of  about  0.37,  which  is  therefore  the  starting 
point  for  the  occurrence  of  capillary  condensation.  It  is  inter¬ 
esting  to  note  that  this  is  in  quite  good  agreement  with  what  is 
displayed  in  Fig.  2,  showing  that  divergency  between  the  adsorp¬ 
tion  and  desorption  isotherms,  which  is  again  the  starting  point 
for  capillary  condensation,  occurs  approximately  for  the  same 
RH. 

Thus,  in  the  model,  for  RH  <  0.37,  no  condensation  occurs 
in  the  pores  of  the  membrane,  and  only  water  adsorption  on 
to  the  surface  of  the  pore  walls  takes  place.  Thus,  proton  con¬ 
ductivity  occurs  only  in  the  thin  water  layers  adsorbed  on  to 
the  pore  walls,  and  its  value  is  relatively  small  (see  the  curve 
representing  the  conductance  g  obtained  at  RH  =  0.2  reported 
in  Fig.  10).  On  the  contrary,  when  RH  reaches  0.37,  capil¬ 
lary  condensation  starts  occurring  in  the  smallest  membrane 
pores,  with  radius  of  about  1.05  nm.  When  RH  is  increased 
above  0.37,  the  curve  representing  the  proton  conductance  g 
presents  two  regions,  which  are  clearly  visible  in  Fig.  10  for 
the  curve  obtained  at  RH  =  0.75:  in  the  first  region,  for  small 
pore  radii  (0<r<2.7nm)  the  pores  are  completely  filled  with 
liquid  water,  and  their  conductance  increases  quadratically  as 
a  function  of  the  pore  radius  according  to  Eq.  (10),  reach¬ 
ing  the  maximum  value  of  1.3  x  10-7  £1  for  a  pore  radius  of 
2.7  nm;  for  pore  radii  larger  than  2.7  nm,  which  is  the  criti¬ 
cal  radius  for  capillary  condensation  at  RH  =  0.75,  there  is  a 
sharp  decrease  in  conductance,  due  to  the  fact  that  the  pores 
are  no  more  filled  with  liquid  water.  From  2.7  nm  onwards, 
the  conductance  increases  linearly  with  the  pore  radius  (Eq. 
(11)),  and  for  large  pore  diameters  (around  10  nm)  the  con¬ 
ductance  is  quite  remarkable  again,  since  these  pores,  even 
if  still  empty  in  their  middle,  have  a  significant  amount  of 
water  adsorbed  on  to  their  wall  surface.  In  Fig.  10,  this  type 
of  behaviour  is  visible  for  the  curves  representing  conduc¬ 
tance  versus  pore  radius  obtained  at  RH  =  0.75,  but  also  for 
RH  =  0.6  (in  this  case  the  critical  radius  for  capillary  conden¬ 
sation  is  about  1 .4  nm),  while  for  RH  =  0.95  the  critical  radius 
is  about  1 1.6  nm,  and  therefore  all  the  pores  in  the  range  repre¬ 
sented  in  the  figure  are  completely  filled  with  water,  and  the 
sharp  transition  corresponding  to  rc  is  not  displayed  in  the 
figure. 

Correspondingly,  the  plots  of  i (r(g)  reported  in  Fig.  1 1  show 
that  the  probability  density  is  concentrated  for  very  low  val¬ 
ues  of  conductance  for  RH  =  0.2,  where  it  displays  a  bimodal 
shape,  very  similar  to  the  shape  of  the  pore  size  distribution. 
xfr(jg)  shifts  towards  larger  values  of  conductance  by  increasing 
RH;  indeed,  for  RH  =  0.5,  iff(g)  shows  a  primary  peak  for  val¬ 
ues  of  conductance  of  about  1.15  x  10-8  £2,  corresponding  to 
the  most  probable  pore  sizes  (1.05  nm)  getting  completely  filled 
up  with  liquid  water.  The  same  peak  also  remains  unchanged  at 
RH  =  0.75  and  RH  =  0.95.  All  the  distributions  also  show  a  sec¬ 
ondary  peak;  this  secondary  peak  represents  the  conductance 
of  the  secondary  peak  of  the  pore  size  distribution,  and,  due 
to  its  gradual  filling  up  with  water  due  to  capillary  condensa¬ 
tion,  it  gradually  moves  towards  higher  values  of  conductance 


Fig.  12.  Simulation  ofNafion®  1 17  at  20  °C.  Comparison  between  model  results 
(lines)  and  literature  experimental  data  (symbols),  (a)  Adsorption  isotherms 
(experimental  data  as  in  Fig.  2);  (b)  proton  conductivity  (experimental  data 
from  [17]). 


as  RH  is  increased,  as  represented  in  the  curves  for  RH  =  0.75 
and  RH  =  0.95  of  Fig.  11. 

Fig.  12a  displays  the  results  obtained  from  the  present  model 
in  terms  of  water  adsorption  for  a  Nation®  117  membrane  at 
20  °C  (the  water  desorption  branch  of  the  isotherm  has  not 
been  simulated  since,  due  to  the  occurrence  of  the  Schroeder 
paradox,  two  different  values  of  water  uptake  are  obtained  exper¬ 
imentally  for  RH=  1,  which  cannot  be  obtained  by  the  model). 
The  agreement  between  simulation  data  and  experimental  liter¬ 
ature  reported  in  Fig.  12a  is  good.  In  parallel,  the  conductivity 
of  the  membrane  has  also  been  calculated,  and  the  results  are 
reported  in  Fig.  12b  together  with  experimental  literature  data; 
the  agreement  is  also  good  in  this  case. 

Fig.  12  shows  how  the  membrane  water  uptake  and  the 
membrane  conductivity  increase  when  increasing  the  RH.  It  is 
interesting  to  note  that,  even  if  for  RH  =  0.75  most  of  the  pores  of 
the  membrane  are  completely  filled  with  water,  no  ‘level  off’  sit¬ 
uation  is  reached  for  the  average  membrane  conductivity  (please 
note  that  the  y-scale  of  Fig.  12b  is  logarithmic).  This  is  due  to  the 
fact  that,  even  if  at  RH  above  0.75  there  is  a  very  small  number 
of  pores  which  are  still  empty,  nevertheless,  those  pores  are  the 
largest  ones,  and  thus,  when  filled  with  water,  they  give  a  sig¬ 
nificant  increment  to  the  overall  conductivity.  This  is  displayed 
by  Fig.  12,  which  shows  a  continuous  increase  in  conductivity 
up  to  RH  close  to  1 . 
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The  model  can  also  be  applied  to  higher  operating  temper¬ 
atures  (in  any  case,  below  100  °C  at  the  operating  pressure 
of  1  atm);  since  the  pore  size  distribution  of  Nation®  mem¬ 
branes  changes  with  temperature,  a  proper  pore  size  distribution 
needs  to  be  known  for  the  operating  temperature  under  conside¬ 
ration. 

4.  Conclusions 

In  this  paper,  a  model  is  proposed  for  the  simulation  of 
Nafion®  proton  conductivity,  where  it  is  assumed  that,  when 
starting  from  a  dry  membrane  and  by  increasing  RH,  membrane 
pore-filling  is  due  to  water  adsorption  on  the  pore  walls  coupled 
to  capillary  condensation;  for  each  value  of  RH,  some  of  the 
pores  (the  smallest  ones)  are  completely  filled  with  water,  while 
the  largest  pores  contain  only  a  thin  layer  of  water  adsorbed  on 
their  pore  surface.  The  conductance  of  each  pore  is  thus  very 
different  depending  on  its  radius,  and  the  overall  conductivity  of 
the  membrane  depends  on  the  percolation  through  the  pore  net¬ 
work.  In  this  paper,  a  model  of  Nafion®  membranes  is  proposed, 
which,  starting  from  literature  data  for  the  pore  size  distribution 
[6],  evaluates  the  occurrence  of  capillary  condensation  in  the 
pores  of  the  membrane,  calculates  their  proton  conductance,  and 
finally  the  overall  proton  conductivity  of  the  membrane  taking 
into  account  percolation  within  the  pore  network  through  the 
effective  medium  approximation. 

The  simulation  results  of  proton  conductivity  of  a  Nafion® 
117  membrane  at  20  °C  show  good  agreement  with  literature 
experimental  data,  even  when  varying  the  RH  operating  condi¬ 
tions. 

The  model  is  independent  of  membrane  thickness,  and  thus 
the  results  presented  in  this  paper  hold  for  Nafion®  membranes 
of  any  type,  provided  that  their  pore  size  distribution  is  the  same 
as  that  used  in  this  paper  [6],  which  was  obtained  specifically 
from  Nafion®  117  membranes.  The  model  can  also  be  applied 
to  higher  operating  temperatures  (in  any  case,  below  100  °C  at 
the  operating  pressure  of  1  atm);  since  the  pore  size  distribution 
of  Nafion®  membranes  changes  with  temperature,  a  proper  pore 
size  distribution  needs  to  be  known  for  the  operating  temperature 
under  consideration. 

This  type  of  model  could  be  applied  (provided  that  the  pore 
size  distribution  is  known)  also  to  other  proton  conducting  poly¬ 
mers  formed  by  a  backbone  with  sulphonic  groups  (or  other 
protogenic  groups)  linked,  such  as  those  polymers  synthesised 


by  the  Dow  Chemical  Company  and  the  Asahi  Chemical  Com¬ 
pany. 
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